Stability analysis of roadbed under flood scouring

Soil roadbed along the river suffers from water erosion at the bottom and collapse at the top under flood scouring, which leads to the suspension of upper pavement slab. In order to ensure the safety of soil roadbed along the river, this study explored the development mechanism of soil roadbed damage by flood in actual cases, and proposed the evolution process of instability under roadbed scouring. The stability law of roadbed along the river under flood scouring was analyzed, and the stability safety factor was corrected to analyze the sensitivity of water depth, flow rate, river bending angle and stability safety factor K in working conditions. The sensitivity of width and height of soil roadbed after flood scouring to water depth, flow velocity, river bending angle was investigated. Moreover, numerical simulation was carried out to determine the displacement nephogram and maximum shear stress nephogram of soil roadbed along the river under the conditions of road surface and roadbed load, vehicle loading or constant change of water depth. By comparing the above theories and engineering cases, the water damage mechanism of soil roadbed along the river was further verified.

The complex topography restricts the planning and design of mountain roads.Take central and western China for example, research data and field survey results show that there are a large number of roads along the river in the western mountainous areas.With the promotion of the West China Development strategy, there have been huge demands for new roads along the river.Roads along the river in mountainous areas are mostly semifilled and semi-dug roadbed, which is vulnerable to long-term scouring by rivers and direct impact by floods or mudslides, resulting in disasters such as roadbed damage, pavement suspension or fracture.This seriously affects the traffic operation of roads in mountainous areas and poses great security risks to local economy and people's lives and property.At present, the flood damage of highway along the river has not been cured, instead, it has become more and more frequent.Therefore, it is of great significance to study the evolution law of water erosion damage of roadbed along the river and implement immediate protection.
Among the studies of roadbed instability due to river water erosion, Chen 1 solved the formula of critical water erosion groove radius of roadbed gap according to the limit equilibrium theory to determine the gap formation mechanism.Zhao et al. 2 carried out numerical simulation analysis to simulate the change of roadbed along river under flood scouring.According to Liang Dan et al. 3 , the erosion and waterout mechanism of roadbed along the river is essentially the interaction process between the erosion resistance of roadbed rock and soil and the impact force of roadbed near the wall.Kawajiri et al. 4 explored the ground erosion and roadbed changes supporting piers by carrying out large-scale open channel model tests.Kurdistani et al. 5 studied the effects of structural geometry, riverbed material and river hydraulic conditions on the downstream scour pattern of stepped weir.Campbell et al. 6 investigated the integrated geophysical technology for underwater detection of bridge foundation related scour and erodible scour filler.Daneshfaraz et al. 7 conducted numerical simulation of river scouring pile groups.Choudhary et al. 8 developed ANFIS-and GEP-based model for prediction of scour depth around bridge pier.
Based on the above analysis, it can be seen that the studies of erosion of roadbed along the river in many countries focus on soil slope erosion test (Khan et al. 9 , Guan et al. 10 ), numerical simulation (Nguyen et al. 11 ) and result analysis (Sharafati et al. 12 ).Through the evolution of water erosion groove formed by flood scouring roadbed toe, this study analyzed the concrete mechanical instability evolution process of roadbed scouring in operation stage, and explored the sensitivity law of roadbed stability safety factor and related factors.It is an urgent

Scouring characteristics of roadbed along the river
The studies of scouring mechanism of soil roadbed along the river focus on the following four aspects: the evolution mechanism of soil roadbed scouring instability, the main controlling factors of soil roadbed scouring, the easily damaged sections of soil roadbed scouring and the disaster time.First, the scouring process of soil roadbed along the river begins with sediment incipient motion.It is considered that the sediment particles or slope soil particles at the slope toe are washed away by river after sediment incipient motion.The porosity of the slope soil increases, the soil roadbed infiltrates, and the slope strength index decreases, which result in the instability and collapse of soil roadbed slope, the scouring damage of soil roadbed and the suspension or collapse of pavement slab.Second, there are three main factors influencing soil roadbed scouring along the river: soil roadbed conditions, water flow conditions and watershed boundary conditions.Third, the most unfavorable scouring position of soil roadbed along the river is at the concave bank of river bend.Fourth, the damage of soil roadbed along the river mainly occurs at the flood scouring stage at the soil roadbed side.
As shown in Fig. 1a, flood scouring leads to the lack of soil roadbed at the concave bank and obvious suspension of pavement slab.Figure 1b is a plan view of roadbed washout at the concave bank, in which Sect. 1-1 shows four stages of soil roadbed lack under flood scouring, as shown in Fig. 2.
(1) Flood scouring stage: In flood season, the velocity and discharge of water flow increase rapidly, and there are a large amount of sediment and gravel in flood, which repeatedly washes away the soil roadbed slope toe.After the scouring action of the main flood circulation 123 reaches the critical state of sediment incipient motion, the soil at the soil roadbed slope toe will be eroded by scouring water, as shown in Fig. 2a.
(2) Formation and expansion stage of water erosion groove: Under normal river conditions, the strength of soil roadbed will gradually decrease after long-term immersion of the soil roadbed slope toe, and there will be potentially dangerous sliding surfaces in soil roadbed soil.Assuming that the rising flood in the flood season constantly washes away the toe of the water-eroded soil roadbed, the initial arc-shaped water-eroded groove a-b-c is locally derived from the soil roadbed.Four potential sliding surfaces I, II, III and IV are developed inside the roadbed, as shown in Fig. 2b.
(3) Soil roadbed sliding instability stage: After the soil roadbed is eroded by flood, the soil at the top of the water erosion groove will expand from a-b-c to def with the rise of flood level under the action of water flow and gravity.In the meanwhile, the upper soil will slide along the sliding surface I. Water erosion groove expands to h-i-g under continous flood scouring, and the upper soil slides along the sliding surface II, as shown in Fig. 2c.
(4) Suspension and instability stage of pavement slab: the water erosion groove continuously expands to the vicinity of the limit state k-m-n.With the increase of the suspension ratio, when the pavement slab reaches the critical state of fracture and failure, fracture and failure will occur, as shown in Fig. 2d.

Stress analysis of soil erosion at the roadbed slope toe
Soil roadbed along the river is often washed by the river, and the sliding force on the slope surface caused by river impact and gravity at the soil roadbed slope toe is related to soil roadbed slope.The stress at the soil roadbed slope toe in Fig. 1a is analyzed, as shown in Fig. 3.

Load analysis of soil at the top of roadbed slope toe
Assume that the inclination angle of roadbed slope is β,the included angle between water flow and the horizontal axis of roadbed slope is δ, the impact drag force is F D ′, the lifting force is F L , the weight of soil is G′, the vehicle load is p i , the weight of pavement slab is q i , and the viscous force of sediment is C k , as shown in Fig. 3a,b.The resultant force acting on the slope surface together is as below: www.nature.com/scientificreports/ The corresponding anti-sliding force is: The starting conditions of sediment at the roadbed slope toe are as follows According to the above formula, the impact drag force of roadbed slope toe is as below: Analysis of soil load on the bottom surface of roadbed slope toe For the bottom surface of soil roadbed slope toe, as shown in Fig. 3c,d, the corresponding sliding resistance is as follows: The impact drag force of soil roadbed slope plane is F D , and the condition of sediment incipient motion is: Set τ c as the incipience drag force of soil roadbed toe bottom surface and τ c ′ as the incipience drag force of roadbed toe top surface.
Substitute formulas (4)-( 6) into formula (7)   When δ = 0, that is, the water flow is orthogonal to the inclined direction of soil roadbed slope toe, the sediment incipience conditions of soil roadbed slope along the river are as follows: F L = 0 , and the lifting force is ignored: According to Formulas ( 6)- (10), under the same water flow, the shear force required for sediment incipience on the soil roadbed slope is less than that on the bottom surface and sediment incipience is more likely to occur on the top surface of the slope toe.Therefore, the river water erosion drives the top surface of the soil roadbed slope toe earlier and easier compared with the bottom surface.

Stress analysis of soil falling on the upper part of roadbed water erosion groove
After the water erosion groove at the bottom of soil roadbed soil is missing, the stress on the upper soil changes, and the soil gradually falls.In videos of flood scouring the soil roadbed along the river, the soil on the upper part of the soil roadbed is irregular (the dotted lines in Fig. 2b-d), falls vertically to the bottom and is washed away by the flood, and the rest is deposited.
Figure 4a,b analyze the stress.Assume that the horizontal dip angle of the upper soil roadbed soil is α′, the weight of the soil is G i ′, the vehicle load is p i , the weight of the pavement slab is q i , the shear resistance of the soil is N e , the vertical dip angle is φ, the width of the upper part is B e , and the height is H e : Stress analysis is carried out for the outer and inner soil at the top of roadbed: The main reason for the falling of flaky soil at the top of roadbed is that the bottom part of is eroded by water, softened, slumped and washed away, which will lose its supporting function to the upper flaky soil.When the (1) water erosion and scouring depth of the bottom water erosion groove reaches a certain size before the upper soil falls every time, the upper flaky soil will fall after overcoming the shear resistance between soil particles under the action of gravity.In this process, the depth dimension of the bottom water erosion groove is equivalent to the width Be′ dimension of the upper flaky soil.In addition, due to the environmental factors such as geography, topography, rainfall and seepage, seepage force will be formed, thus affecting the upper soil.Therefore, the flake size of the falling flaky soil is related to the density, material properties, water erosion depth, scouring depth, flake height of the soil and the hanging width, and most of the falling parts are flake structures.
To sum up, from the stress analysis in Sects.3.1 and 3.2, it can be found that the falling position of the roadbed soil along the river caused by water erosion is mainly determined by the water erosion groove of the bottom soil.Therefore, the change of the water erosion groove angle of the soil roadbed will affect not only the integrity and missing direction of the whole soil roadbed, but also the proportion and position of the upper pavement slab suspension.

Overall stability analysis of roadbed with water erosion groove under the gravity of vehicle and pavement slab
Based on field investigation and relevant literature review, it can be found that the mechanism of water-damaged settlement failure of soil roadbed along the river is similar to that of rainfall-induced landslide, and the sliding surface is usually arc-shaped.When the soil roadbed along the river is washed by flood, a water erosion groove is formed at the slope toe, and the upper soil still meets the principle of limit equilibrium in a stable state.Therefore, based on the circular sliding mechanism, this study used Bishop method to analyze the stability of soil roadbed washout along the river.Bishop 13 assumed that: (1) The acting force of soil between strips is only normal force, and there is no tangential force; (2) The sliding soil conforms to the overall moment balance condition; (3) Each block conforms to the polygonal closing condition of force; (4) The anti-sliding safety factor of the bottom sliding surface of roadbed soil strip is the same, which is equal to the average safety factor of the whole roadbed sliding surface.The soil slope per unit length is calculated as a plane problem.Set the landslide surface of roadbed as arc AC, the center of the circle as O and the radius as R. The landslide ABC is divided into several soil strips, and the stress of which is analyzed, as shown in Fig. 5.The stress on the ith soil strip is as below: (1) Weight of soil strip, Gi = γb i h i ; Where, b i and h i are the width and average height of soil strip i; When there is a water erosion groove, Gi = γb r h r ; Where, b r and h r are the residue width and residue height of soil strip i. (2) Shear resistance T fi , effective normal reaction N i ' and pore water pressure u i l i acting on the bottom surface of soil strip i, where, u i and l i are pore water pressure at the midpoint of the bottom surface of soil strip i and the length of sliding arc respectively; The Angle between the center line and the vertical of the soil strip i is α i .(3) Normal force E i and E i+1 and tangential force X i and X i+1 acting on both sides of soil strip i, ΔX i = (X i+1 − X i ).(4) Vehicle weight p i and pavement slab weight q i acting on soil strip i; Moreover, the action points of G i , T fi , N i ′ and u i l i are all at the midpoint of the bottom surface of the soil strip i. www.nature.com/scientificreports/For the vertical direction of soil strip i, according to the balance of force, the following formula can be obtained: When the soil slope has not collapsed and remains stable, the shear strength T fi on the sliding surface of soil strip i is expressed by effective stress, as below: where, c′-Effective cohesion of soil, φ′-Effective internal friction angle of soil, K-Safety factor.
Substitute (12), and N i ′ is where The moment balance analysis of the whole sliding soil body to the center O of the circle is carried out.Since the side wall action moments between adjacent soil strips cancel each other, and the action lines of N i ' and u i l i of each soil strip i pass through the center O of the circle, the following formula can be obtained: Substitute formulas ( 14) and ( 16) into formula (13), and x i = Rsinα, b = b i = l i cosα i .According to Bishop method, the safety factor of roadbed stability under the action of pavement slab and vehicle is as below: ΔX i is estimated to determine K through the successive approximation method.The trial value of X i and E i should meet the equilibrium conditions of each soil strip.∑ΔX i and ∑ΔE i of the whole sliding soil strip should be equal to zero.Bishop proved that the error is only 1% when ΔX i = 0, and formula (17) can be simplified as follows:

Stability analysis of roadbed under straight flood scouring
Considering that the scouring force is parallel to the longitudinal section of the roadbed when the flood scours straightly along the river, as shown in Fig. 6, the moment to the center O of the circle is 0. Therefore, formula (18)  can be used to calculate for calculating the stability safety factor of soil roadbed along the river in the straight section.
Where β is the toe of the slope and β 1 is ∠OCB, as shown in Fig. 8.In this study, in order to calculate the most dangerous state, let F act on the most dangerous position, that is, the toe of soil roadbed slope, and the moment of the center O of the circle can be simplified as M = FR.www.nature.com/scientificreports/According to the momentum theorem Ft = 2mv ⊥ , where v ⊥ is the velocity of water flow perpendicular to the river bank in the concave part of the river, v ⊥ = vcos(θ/2), where v is the river flow velocity, θ is the angle of the river bend, m = ρv = ρQt, and Q is the flow of the river bend, F = 2ρQv ⊥ is obtained, which is substituted into (18), and the stability safety factor at the bend along the concave part of the river road is obtained as follows: When the soil slope is partially flooded.The hydrostatic pressure P 1 acting on the sliding surface below the water level line, the water pressure P 2 on the slope surface, the gravity of pore water and the reaction force G w of soil particle buoyancy are balanced in the hydrostatic state.The action line of P 1 passes through the center O of the circle.According to the moment balance, the moment of P 2 to the center O of the circle and the moment of G w to the center O of the circle cancel each other, as shown in Fig. 8. Therefore, under the condition of hydrostatic state, the influence of water pressure on roadbed surface on the sliding soil can be calculated according to the buoyancy of the sliding soil below the still water surface, that is, the weight of the underwater soil strip is calculated according to the effective weight.The calculation formula of the stability safety factor is the same as formula (19).
When it rains continuously, the mountain surface water flows downwards and seeps into the river through the roadbed.The water level in the slope is higher than that outside the slope, and the water in the slope generates seepage force outward, pointing to the slope surface, as shown in Fig. 8. Assuming that the saturation line is e-f-g, the area of the part (f-g-c) of the sliding soil under the saturation line is A w , and the total seepage force D on the soil is as below: where, J-Seepage force on unit volume of soil (KN/m 3 ); i-The average value of hydraulic gradient in the area A w under the saturation line, i is assumed to be equal to the gradient of the connection line of the saturation line f-g.γ w -Saturated water severity.
The action line of the seepage resultant force D is located at the centroid of area f-g-c, assuming that the action direction is parallel to f-g, and the arm of force of D to the center O of the circle is recorded as r.Therefore, considering the seepage force, formula ( 19) is substituted and the formula for analyzing the stability safety factor of roadbed along the river by Bishop method is as follows: Homogeneous, laminar and stable gradual flow are assumed according to the hydraulics [15][16][17] , as shown in Fig. 9.
The upper water level is H 1 , which can be approximated as the height H of roadbed along the river in this study.The height of lower water level is the height H 2 of river water level, and the horizontal projection length of roadbed bottom is L. According to Darcy's law, the average velocity of roadbed across the river is as below: The seepage flow is q s : The variables are separated from the above formula, and the upstream surface ( x = 0, z = H 1 ) H is integrated to the downstream surface ( x = L, z = H 2 ): According to the saturation line formula (28), Aw can be obtained.

Prediction and analysis of roadbed erosion by river
In order to calculate the stability safety factor of the residue roadbed soil of the highway along the river in water erosion, the self-weight of soil strips is recorded as G ir = γb ir h ir , where, b ir and h ir are the width and height of soil strip i of the residue roadbed soil after river water erosion, respectively, and they are substituted into Formula (30).
In the state of smooth revetment, as shown in Fig. 10, the roadbed along the river is partially scoured by flood.Based on the scour depth calculation formula of smooth revetment in Code for River Regulation Design (GB50707-2011) 14 , the scour depth is as below: According to the stable sliding surface obtained by the Bishop method, combined with formula (32), the residue stable roadbed width is obtained as follows: Based on formula (33), the calculation formula of h r corresponding to the critical soil strip for b r is as follows: According to the vertical distance h d between the roadbed drawn by the Bishop method and the vertex of the center line of unstable soil strip i, it is concluded that the height of the residual stable roadbed soil after river scouring and water erosion is H r .
The height of the roadbed bottom eroded by river water is as follows: In order to ensure the safety of highway operation and protect the stability of pavement slab and roadbed, the limit value of passing vehicles should be obtained 18 .Therefore, according to formula (36), the vehicle load P l is as below:

Case analysis
The bend near K4 + 210 third-class highway in Shizhu County, Chongqing is taken as an example.This section is a semi-filled and semi-dug soil roadbed, with soil roadbed height H = 5 m and pavement width b = 4.0 m.The section along the river is a homogeneous soil-rock mixture slope.For the semi-filled and semi-dug roadbed, the medium-dense pebble soil, gravel soil, rounded gravel soil and breccia soil are used for filling, and the particle size of sand is less than 150 mm.The slope is a Class IV strong wind fossil slope, which is affected by groundwater 19,20 .The filling slope is 1:1.According to the Bishop's method, soil strips are divided into 15 strips, and the stability safety factor K of soil slope is calculated by Formulas (18-37).The changing height and width of scouring water erosion at the bottom of roadbed along the river and the load limit value of vehicles passing on the road are calculated and predicted.Physical and mechanical indexes of roadbed filler are shown in Table 1.

Sensitivity analysis of water depth, river bending angle and stability safety factor
In the depth of the water, river bending angle and the stability safety factor of sensitivity analysis, working condition of parameter selection are as follows: Roadbed slope along the river, the slope toe slope ratio is 1:1, the flood depth is H 2 = 1.6 m, 2.0 m, 2.4 m, 2.8 m, 3.6 m and 4.0 m; respectively, and the water flow velocity is 1.8 m/s, 2.0 m/s, 2.2 m/s, 2.5 m/s, 2.8 m/s and 3.1 m/s.The bending angle at the concave bank of the smooth revetment river is θ = 90°,105°,120°,135°,150°,165°,180°; as shown in Figs.11, 12, 13.

Sensitivity analysis of flow velocity, river bending angle and stability safety factor
In the flow velocity, river bending angle and stability safety factor of sensitivity analysis, working condition of parameter selection are as follows: Roadbed slope along the river, the slope toe slope ratio is 1:1, the flood depth is 1.Relevant parameters of roadbed and river are selected for the above working conditions.After the design and construction of roadbed is completed, although the design K value has been determined, the sensitivity analysis of the demand value of stability safety factor for the above common working conditions of roadbed along the river under flood and water erosion can be taken as a reference to guide the safety and protection of roadbed according to the following analysis conclusions:    www.nature.com/scientificreports/ (1) When the flood level H 2 changes from 1.6 m to 4.0 m, the demand value of stability safety factor K of soil roadbed increases with the increase of the water depth; (2) When the river bend angle θ changes from 90° to 180°, the demand value of stable safety factor K decreases with the increase of the river bend angle; 3) When the water depth and velocity are fixed, the required value of stability safety factor K decreases with the increase of the river bending angle.(4) When the water depth is fixed, the required value of stability safety factor K increases with the increase of river velocity, and the increase rate decreases with the increase of the river bending angle.Particularly, the required value of stability safety factor K is constant for straight river section.( 5) When the flow rate is fixed and the river bend angle is the same, the demand value of stability safety factor K increases with the increase of the water depth; (6) When the flow rate is fixed and the bend angle of the river becomes larger, the required value of the stability safety factor K decreases with the increase of the water depth.( 7) By comparing the required value K with the original design value, it can be judged whether the soil roadbed along the river exceeds the original design stable state under the corresponding working conditions.When the required value of K is smaller than the original design value, the soil roadbed along the river is safe under the corresponding working conditions.When the required value of K is close to the original design value, it is in a critical state, so it is suggested to take corresponding reinforcement measures; As the required value of K is larger than or equal to the original design value, the soil roadbed along the river is likely to be damaged, and corresponding preventive measures should be taken.Therefore, by comparing the required value of K with the original design value, it can be timely and accurately judged whether the soil roadbed along the river is in a safe state under the corresponding working conditions to take corresponding reinforcement and damage prevention measures, so that pre-judgment and decision-making for disaster prevention and mitigation can be realized in advance.In this working condition, the common parameters of soil roadbed and river are selected.When the soil roadbed along the river at different river bending angles is scoured under different river flow rates, the width and height of residue roadbed soil change.Through the sensitivity analysis of flow rate, river bending angle and river water level on the width and height dimensions of residue soil roadbed along the river, the following conclusions can be drawn: (1) When θ ≤ 120°, the erosion damage of the slope and toe of soil roadbed along the river increases with the increase of water depth and velocity, and the corresponding values b r and h r decrease.(2) When θ = 90°, the slope and toe of soil roadbed along the river suffer from the strongest scouring destructive power, the corresponding b r and h r decrease the most rapidly, and their corresponding scouring change rate are also the largest.However, with the increase of bending angle, the change rate of the scouring decreases.(3) The river bending angle is the most obvious factor that affects the slope and toe of soil roadbed along the river.Particularly when the river bending angle is 90°, slope protection measures should be taken.(4) When the river bending angle is θ = 135°, it can be seen from the calculation results that v = 2.2 m/s is the scouring boundary value of soil roadbed slope and toe along the river.when v ≤ 2.2 m/s, b r and h r increase

Numerical simulation analysis of roadbed along the river under flood scouring
The model was established by FLAC3D: The size of the pavement slab and the soil roadbed is 5 m × 4.0 m × 0.3 m and 5 m × 4.0 m × 3 m respectively, and the slope is 1:1, as shown in Fig. 22.The material parameters of this model are listed in Table 2.The model adopts automatic grid division, with 131,327 units and 117,677 nodes.The scouring simulation of soil roadbed slope along the river was carried out in turn under the river water depth of 1.6 m, 2.0 m, 2.4 m, 2.8 m, 3.2 m, 3.6 m and 4.0 m, the water flow velocity of 2.2 m/s, and the river bend angle of 90°.The changing trend of displacement and stress of roadbed soil was verified when the road vehicle was driven by a single wheel with 100KN against the outside of roadbed slope and the pavement slab.Through numerical simulation of the displacement nephogram of soil roadbed along the river under the condition of only road surface and soil roadbed load, vehicle loading and continuous water depth changes, as shown in Fig. 23, the following change law can be found: (1) Under road surface and soil roadbed load, the displacement nephogram shows that the soil roadbed soil under the road surface is subjected to its gravity, resulting in vertical settlement displacement.With the increase of the soil roadbed depth, the settlement displacement decreases.This is because the soil at the bottom is compacted and the soil at the top is vertically settled.(2) When the vehicle is driving on the road surface, the settlement displacement of roadbed soil under the road surface increases.With the increase of roadbed soil depth, the settlement displacement is higher than that when there is only road surface and roadbed load.(3) When flood scours roadbed soil, the scour displacement at the soil roadbed slope toe changes and water erosion grooves appear.(4) The water erosion groove at the soil roadbed slope toe increases with the increase of flood water level, and the displacement of roadbed soil slope also increases with the increase of water level, which is connected with the soil at the bottom of pavement slab and tends to landslide in the form of sliding surface.
Through numerical simulation of the maximum shear stress nephograms of soil roadbed along the river under the conditions of only road surface and roadbed load, adding vehicle loading and continuous water depth change, as shown in Fig. 24, the following changing law can be found: (1) When only the road surface and soil roadbed are loaded, the maximum shear stress nephogram shows that the shear stress on the bottom roadbed soil is the largest, and it gradually decreases from bottom to top.This is because the bottom roadbed soil is compacted under the action of gravity.(2) When there are vehicles driving on the road, the shear stress of the road surface is the largest, followed by the overall shear stress of the roadbed soil, which is large and uniform, and the shear stress of the upper part of the roadbed slope is small.(3) When the flood scours roadbed soil, the maximum shear stress at the roadbed slope toe changes, and the water erosion grooves appear.(4) Under the action of maximum shear stress, the stress nephogram of soil roadbed slope changes, there is erosion and loss of roadbed soil.With the increase of water level, the height of loss gradually increases, and the maximum shear stress of roadbed soil changes as a sliding surface.
The displacement and maximum shear stress analysis of the above numerical simulation are consistent with the previous stability analysis conclusion.

Figure 1 .
Figure 1.Erosion of roadbed along the river.(a) Water damaged roadbed along the line, (b) River plan.

Figure 2 .
Figure 2. Evolution process of suspended pavement slab.(a) Slope toe of roadbed along the river eroded by concave bank, (b) Water erosion groove of roadbed along the river, (c) Roadbed soil slides along the dangerous sliding surface, (d) Suspension of pavement slab caused by the sliding instability of roadbed soil.

Figure 3 .
Figure 3. Schematic diagram of stress analysis of sediment particles on roadbed slope.(a) Schematic diagram of the top surface of slope toe, (b) Schematic diagram of stress on the top surface of slope toe, (c) Schematic diagram of the bottom surface of slope toe, (d) Schematic diagram of stress on the bottom surface of slope toe.

Figure 4 .
Figure 4. Schematic diagram of stress analysis of soil at the top of roadbed.(a) Stress analysis of outer soil falling at the top of roadbed, (b) Stress analysis of inner soil falling at the top of roadbed.

Figure 5 .
Figure 5. Diagram of Bishop calculation of roadbed with water erosion groove under vehicle loading.

Figure 6 .
Figure 6.Diagram of straight flood scouring roadbed along the river.

Figure 7 .
Figure 7. Diagram of bend flood scouring roadbed along the river.

Figure 8 .
Figure 8. Calculation of roadbed stability during water seepage.

Figure 11 .
Figure 11.Sensitivity analysis of different H 2 and θ to K when V = 1.8 m/s, 2.0 m/s.

Figure 12 .Figure 13 .
Figure 12.Sensitivity analysis of different H 2 and θ to K when V = 2.2 m/s, 2.5 m/s.

Figure 14 .
Figure 14.Sensitivity analysis of different v and θ to K when H 2 = 1.6 m/s, 2.0 m/s.

Figure 15 .
Figure 15.Sensitivity analysis of different v and θ to K when H 2 = 2.4 m/s, 2.8 m/s.

Figure 16 .
Figure 16.Sensitivity analysis of different v and θ to K when H 2 = 3.2 m/s,3.6 m/s.

Figure 17 .
Figure 17.Sensitivity analysis of different v and θ to K when H 2 = 4.0 m/s.

Figure 18 .
Figure 18.Sensitivity analysis of different v, H 2 and b r of the roadbed along the river when θ = 90°, 105°.

Figure 19 .
Figure 19.Sensitivity analysis of different v, H 2 and b r of the roadbed along the river when θ = 120°, 135°.

Figure 20 .
Figure 20.Sensitivity analysis of different v, H 2 and h r of the roadbed along the river when θ = 90°, 105°.

Figure 21 .
Figure 21.Sensitivity analysis of different v, H 2 and h r of the roadbed along the river when θ = 120°, 135°.

Figure 22 .
Figure 22.Diagram of numerical simulation model.

Figure 23 .
Figure 23.Displacement nephogram of roadbed pavement.(a) Displacement cloud image under the action of gravity on subgrade road surface, (b) Displacement cloud image under the action of vehicle and gravity, (c) Displacement cloud image with water erosion water groove at water depth of 1.6 m, (d) Displacement cloud image with erosion groove at water depth of 2.0 m, (e) Displacement cloud image with water erosion, groove at water depth of 2.4 m, (f) Displacement cloud image with water erosion groove at water depth of 2.8 m, (g) Displacement cloud image with water erosion at water depth of 3.2 m (h) Displacement cloud image with water groove erosion groove at water depth of 3.6 m, (i) Displacement cloud image with water erosion groove at water depth of 4.0 m.

Table 1 .
Physical and mechanical indexes of roadbed filler.

Table 2 .
Physical and mechanical indexes.